Multiphoton Ionization as a "Clock" to Reveal Molecular Dynamics with Intense 

Short X-ray Free Electron Laser Pulses 
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We investigate molecular dynamics of multiple ionization in N2 through multiple core-level pho- 
toabsorption and subsequent Auger decay processes induced by intense, short X-ray free electron 
laser pulses. The timing dynamics of the photoabsorption and dissociation processes is mapped 
onto the kinetic energy of the fragments. Measurements of the latter allow us to map out the 
average internuclear separation for every molecular photoionization sequence step and obtain the 
average time interval between the photoabsorption events. Using multiphoton ionization as a tool 
of multiple-pulse pump-probe scheme, we demonstrate the modification of the ionization dynamics 
as we vary the x-ray laser pulse duration. 

PACS numbers: 32.30. Rj, 82.80.Ej, 33.20.Rm, 33.60.+q 



Substantial efforts have been made both experimen- 
tally and theoretically to explore molecular dynamics 
during and after the interaction of the molecule with ra- 
diation, i.e., probe the electrons and nuclei behavior sub- 
sequent to absorption of photons [l]-[i| . Among the num- 
ber of factors determining this dynamics, the wavelength 
and intensity of the radiation are particularly important. 
Recent development of free electron lasers (FELs) has 
provided a novel source of short X-ray pulses with much 
higher intensity than available before |5(. However, first 
experiments on molecular dynamics in this, hitherto un- 
available range of parameters, have revealed that our 
understanding of the involved physics and chemistry is 
largely qualitative. We are lacking accurate models for 
quantitatively predicting the behavior of molecular sys- 
tems that are subjected to intense X-rays interacting di- 
rectly with atomic inner shells. The motivation of the 
present study is to remedy this situation by investigat- 
ing molecular dynamics and fragmentation, subsequent 
to core ionization and core-hole relaxation, and by pro- 
viding a model and quantitative description of the in- 
teraction. These studies are of fundamental importance 
to molecular physics, physical chemistry and solid state 
physics. Furthermore, we expect that, with the rapid de- 
velopment of the FEL sources around the world, the need 
and the ability to accurately model the early phase of the 
molecular target breakup (that inevitably follows the in- 
tense X-ray pulse) becomes an essential part of many 
studies. 

The interaction of the intense X-ray FEL pulse with 
a molecule starts in the same way as for a lower inten- 



sity X-rays, i.e., with single-photon absorption and the 
ejection of a core electron. In lighter elements such as 
nitrogen the resulting core hole is filled by Auger decay 
in a few femtoseconds (6.4 fs for N2), with the ejection of 
a second electron (Auger electron). We will refer to this 
sequence as the photoionization- Auger (PA) cycle. Using 
conventional X-ray sources, the likelihood of the target 
molecule absorbing another photon before the dissocia- 
tion is complete is negligible. In contrast, the photon flux 
in the FEL pulse is so high that the dissociating molecule 
has a chance to absorb several X-ray photons in sequence, 
so that multiple PA cycles (a PAPA. . . sequence) can take 
place in the molecular system within just one X-ray FEL 

puisc mnj. 

In this work, we present an alternative to conventional 
pump-probe schemes @, EI"13|, taking advantage of the 
femtosecond timescale of the multiple PA cycles in in- 
tense FEL pulses. With suitably chosen pulse inten- 
sity, the average interval between the PA cycles can be 
matched to the dissociation timescale, so that the ob- 
server follows the progress of the molecule along the re- 
pulsive potential energy surfaces (see Fig. [TJ. We can 
regard the multiple PA cycles as an 'intrinsic" novel 
multiple-pulse pump-probe mechanism where each sub- 
sequent cycle probes the development of the dynamics. 
Both the pumping and multiple probing come from a 
single X-ray pulse by sequential absorption. In the de- 
scribed multiphoton scheme, applicable to any ionization 
with intense X-ray-FEL pulses, the pump-probe delays 
are controlled by varying the FEL parameters. 

In contrast to the pump-probe schemes based on two 



2 



external pulses, the timing in the present scheme is nec- 
essarily of statistical nature, since it is determined by 
intrinsic quantum events. At each PA cycle, the kinetic 
energy (KE) of the emitted electrons reflects the molec- 
ular potential curves of the relevant charge states at a 
particular internuclear separation (Fig.Q}. The eventual 
kinetic energy release (KER) to the fragment ions, on 
the other hand, is determined by the entire dissociation 
dynamics. By combining these various observables with 
quantitative modeling, one can extract the averaged tim- 
ing information and reconstruct molecular dynamics. 

Thus, on one hand the PA processes provide a means 
for studying the electronic and structural dynamics in 
the X-ray-molecule interaction. And on the other hand 
it is also an atomic intrinsic tool of quantum control of 
the interaction, since by manipulating the X-ray beam 
temporal profile one can modify the dynamics of sequen- 
tial multiple ionization, opening or closing PA pathways 
to select desired outcomes [1, [l^ . 

In the present study, we examine the measured multi- 
ple ionization sequences of N2 molecules with X-ray FEL 
pulses at equal pulse energies and various pulse dura- 
tions, which result in different dynamics profiles of the 
interaction. We model the scenario of the interaction 
and quantitatively map the time dependence onto the 
average KE of the fragment ions. We extract the aver- 
age time interval between photoabsorption events At aV g 
as well as the corresponding internuclear separation (R) 
at which each step occurs in an ionization sequence. As 
expected, we observe that At avg decreases with increas- 
ing intensity: a) at equal pulse energy which leads to a 
certain final charge state, At avg decreases with decreas- 
ing pulse duration, e.g., At avg decreases from 11-23 fs to 
7-12 fs as the pulse duration decreases from 280 fs to 80 
fs; b) at equal pulse duration, At avg decreases with in- 
creasing pulse energy, e.g., At avg decreases from 12-23 fs 
to 7-11 fs in X-ray spatial regions with increasing pulse 
energies where N 3+ to N 7+ are preferentially produced. 
This work takes advantage of sequential PA processes to 
reveal the dynamics timing by developing a model, which 
is validated by the measurement of the fragment charge 
states. It shows the ability to quantify the correlation be- 
tween ionization rates and KERs, and the possibility to 
resolve individual photoabsorption step in multiphoton 
ionizations. Also, this work presents the modification of 
At avg with manipulation of the X-ray pulse duration and 
shows the variation of the ionization timing dynamics due 
to different pulse intensities within a single pulse. 

In the current experiment, the first PA step produces 
molecular dications N^ 4 " in either the ground state or ex- 
cited states, with vacancies in outer molecular orbitals. 
About 80% of the electronic states of the dications that 
are formed in the Frank-Condon region, are dissocia- 
tive 1J, [15J. Throughout this work, we will refer to 
the dissociation process N^" 1 "™ ->N" + +N m+ as the (n,m) 
channel. While the molecule breaks up, other PA cycles 



can take place locally at each of the atoms and thus the 
dissociation process may jump from the (1,1) to the (3,1) 
and then to the (3,3) channel. For a stable molecular di- 
cation the dissociation can be started in the (2,2) chan- 
nel following the second PA process, then jumps to (4,2) 
and (4,4) channels by subsequent PA's. Similarly, higher 
charge states can be populated with further PA steps. 
The ionization sequences can be primarily defined by the 
final charge of the atomic fragments, but the sequences 
have common intermediate ionization steps. (See sum- 
mary of possible pathways in Ref. 



Starting with the ground state of neutral N 2 , each pho- 
toionization step will lift the newly produced molecular 
ion to a potential curve of a higher charge state, as illus- 
trated schematically in Fig. [TJ Each transition occurs at 
a different time and, once the dissociation starts, at a dif- 
ferent R. For an ionization sequence, At avg is determined 
by the inverse rate of the core ionization, which is given 
by the product of the core ionization cross-section and 
the X-ray flux. Hence, the higher the X-ray-FEL inten- 
sity, the shorter At avg . This gives dissociation less time 
to proceed between ionizations and hence corresponds to 
a smaller R for each PA step. The smaller the R at 
which an ionization occurs, the stronger the influence of 
the neighboring atom on the one that hosts the photoion- 
ization process. 
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FIG. 1: (Color online) Schematic of the PA sequences while 
molecules dissociate. The PA sequence shown here as an ex- 
ample is N 2 ->-(l,l)-K3,l)-K3,3). "P": photoionization; "A": 
Auger decay. 

The experiment was performed in the high field physics 
(HFP) end station of atomic molecular and optical 
physics (AMO) hutch at the Linac Coherent Light Source 
(LCLS) H, TniM- The photon energy used was 1100 
(±15) eV, which is above the K-edge of neutral N and 
all charge states of N2 . The nominal pulse energy values 
from the gas detector monitors are measured upstream of 
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the HFP chamber. In the interaction region the actual 
pulse energy is reduced by 65%-85% due to the beam 
transport loss 0, d, [l7| ■ The X-ray beam is focused by 
a pair of Kirkpatrick-Bacz (KB) mirrors to an oval spa- 
tial profile with a major axis of 2.2 /jm and a minor axis 
of 1.2 fxm d, HI- The pulse durations used were 280 fs, 
110 fs, 80 fs and 7 fs, estimated by measuring the width 
of the electron bunch The data shown in this 

work were recorded by an ion time-of-flight (iTOF) spec- 
trometer, which is oriented perpendicularly to the plane 
defined by the polarization and propagation direction of 
the X-ray beam @, & (H . 

To extract the KE of the ion fragments, we use a 
SIMION simulation to obtain a basis set of simulated 
iTOF spectral profiles corresponding to different charge 
states and KEs 0, 23 1. We adjust the weights of the 



basis for the best fit of the experimental spectra. These 
iTOF basis functions are generated in such a way that 
they automatically take into account all the non-trivial 
spectrometer transmission coefficients for different charge 
states and KE values. The average KEs obtained by this 
method give consistent results when using different basis 
sets and pools of initial values for fitting to the experi- 
mental spectra. 

Knowing all the potential energy curves involved in 
a certain dissociation pathway with concurrent PA se- 
quences, and At avg , one can calculate the average KER 
in a channel (m,n). One also needs to know the time 
profile of the intensity of the FEL pulse; in this work 
we have approximated it by a square function [22[ in 
which case the ionization rate is constant. The calcula- 
tion can be performed also in reverse, obtaining KERs 
from experimentally determined KEs and then extract- 
ing R's and times for each preceding PA step. In this 
calculation we: i) assume the dominant contribution to 
the observed N m+ ions is from the (m,m) channel, so that 
the KER is double that of the measured KE; ii) assume 
that there is one dominant pathway leading to ions of a 
certain charge; iii) assume, in the R range of the dissocia- 
tion of core-ionized molecules, core-hole potential curves 
have the same slopes as the potential curves of valence- 
hole states populated by the corresponding subsequent 
Auger decay process, i.e., assume that the Auger decay 
is instantaneous; iv) set At aV g inversely proportional to 
the cross section of the atomic transitions; v) approxi- 
mate potentials of high charge states with Coulomb po- 
tentials (a screening factor is applied to the (2,2) state). 
The validity of these assumptions and approximations is 
discussed in the supplementary material [161 ]. 

With intense X-ray pulses, we observe highly charged 
nitrogen fragment ions as a result of PA sequences, as 
shown in the main panel of Fig. [2] The spectra in Fig. [2] 
are measured with different pulse durations at equal flu- 
ence, except for the 7 fs pulse duration data, for which 
the pulse energy is lower by 84%. The 7 fs spectrum, 
from which the PA sequences are essentially missing and 



which is dominated by processes of molecular nature, is 
shown as a reference (a case where the least molecular 
dissociation is involved). As seen in the main panel of 
Fig. [21 higher X-ray flux results in higher charge states. 
This may suggest a molecular component to the R depen- 
dent atomic core-electron absorption cross section which 
is changing during molecular dissociation. Also, the dy- 
namics leading to the production of bare N 7+ appears to 
saturate. We observe that the pulse duration has a strong 
effect not only on the charge state branching ratios [§], 
but also on the KE of each ion species. In the inset of 
Fig. [2l a wider structure corresponds to a larger aver- 
age KE. For N + and N 2+ , the peak widths are the same 
at the three different pulse durations, but starting from 
N 3+ the peak width decreases with increasing pulse du- 
ration. To quantify this trend, we obtain the average KE 
for different ion fragments at different pulse durations, 
which is shown in Fig. [3] Since molecular ions or atomic 
fragments are detected after all ionization processes are 
over, the final KEs are the sum of contributions from 
all ionization steps involved. The KER shared by the 
ion fragments is determined by the R where the PA or 
photoionization processes occur: the smaller the R, the 
larger the resultant KER. In other words, higher charges 
created at shorter separations convert more repulsive po- 
tential energy into the kinetic energy of the ions. As 
expected and seen in Fig. [3[ for N + and N 2+ , the KEs 
are independent of the pulse duration, since these ions are 
mainly formed by PA processes in bound molecular states 
and the dissociation is not perturbed by further dissoci- 
ation. For ionic fragments with charges >3+ the KE is 
larger for pulses with shorter duration but with similar 
fluence, i.e., a higher intensity (see data with pulse dura- 
tions >80 fs). The observed increase of ion KE for shorter 
duration pulses supports the reasoning that higher pho- 
ton intensity reduces At avg . It is also seen from Fig. [3] 
that, for a given pulse duration the KE increases with 
the ion charge. Even at high charge states, e.g., for N 7+ , 
we see an increase of KE compared with that of N 6+ , re- 
flecting the long range operation of Coulomb repulsion, 
and showing that the dissociation is not complete after 
N 5+ ionization. 

As described previously, with the measured KE and 
using the above model, we extract At avg and the cor- 
responding R for various ionization sequences. This is 
done independently for different final charge states and 
the results are shown in Fig. [4j In this figure, each 
block represents a pathway leading to the correspond- 
ing charge state; slices (or dots) of different colors within 
each block represent the ionization sequences for three 
different pulse durations. In our model, the time of the 
first ionization event is based on the average ionization 
rate, not the start of the FEL pulse, so the absolute start- 
ing times shown in Fig. [4] are not equal. For a given final 
state, the ionization step occurs earlier with a shorter 
pulse than with a longer pulse. For instance, the R (and 
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FIG. 2: (Color on line) Ion spectra at different pulse dura- 
tions. The pulse energy is 1.6-1.9 mj (nominal values) for 
pulse durations >80 fs. The pulse energy for the 7 fs pulse is 
0.2-0.3 mj (nominal values). Main panel: spectra are normal- 
ized to the total ion signal of all charge states. Insets: spectra 
for different charge states are normalized to the maximum of 
each spectral peak. 
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FIG. 3: (Color on line) Average KE of ion fragments at dif- 
ferent pulse durations, obtained from the spectra in Fig. [2] 

At avg ) to reach N 3 + with 80 fs, 110 fs and 280 fs arc 
1.73±0.09 A(12 fs), 1.90±0.15 A(15 fs) and 2.37±0.20 
A(23 fs). We also find that, for a given pulse duration, 
the time to reach a specific intermediate charge state is 
different for ionization sequences leading to different fi- 
nal charge states. For example, using the KE of N 3+ , 
N 5+ , N 6 + and N 7 + at 280 fs, R (and At avg ) for the tran- 
sition (1,1)->(3,1) arc found to be 2.37±0.20 A(23 fs), 
1.87±0.15 A(14 fs), 1.73±0.09 A(12 fs) and 1.69±0.08 
A(ll fs), respectively. We thus see a decreasing R and 
At avg for a given intermediate transition in ionization se- 



quences associated with increasing final ion charges, as 
they are produced preferentially at increasing energies 
within a single pulse (see inset of Fig. 2]). 

The spatial intensity distribution in the focus leads to 
regions of high and low ionization rates and the mea- 
sured KE is the convolution of these regions. Our ability 
to perform "pulse length resolved" spectroscopy of disso- 
ciating N2 molecules reveals that the involved molecular 
dynamics are influenced by both the spatial and tempo- 
ral profile of the FEL pulse. Furthermore, the accuracy 
of the determined timing and R for ionization steps de- 
pends on the potential energy curves used in the model. 
One could retrieve the information of the potential en- 
ergy curves through measurements of the timing and R 
which could be achieved in a standard pump-probe ex- 
periment. 




FIG. 4: (Color on line) R as a function of time for different 
PA pathways at different pulse durations. Each colored dot 
represents a photoabsorption event. The gray dots represent 
the projection of the event onto the channel-time plane. The 
channel axis labels indicate the states after the first PA se- 
quence and the final states. The potential curve used 
in the calculation is the 1 A g state from Ref. [24|. Inset: 
schematic of the spatial intensity profile of the FEL beam; 
the labels correspond to the letter labels in the main panel, 
indicating that higher charge-state ions are preferentially pro- 
duced in the region with higher pulse intensity. 

In summary, we investigate the ionization dynamics of 
N2 molecules through sequential multiphoton core-level 
ionization and subsequent Auger decay. Qualitatively, 
we expect and observe the molecular dynamics varying 
as a function of the X-ray intensity. With ion modeling, 
we are able to quantify the correlations between pulse 
durations, photoabsorption rates and fragment KEs. We 
obtain the average ion KE versus charge state and extract 
At avg and the corresponding R. With the multiphoton 
ionization scheme, we demonstrate the modification of 
the dynamics by varying the X-ray temporal profile. To 
completely disentangle the photoionization and the fol- 
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lowing Auger decay and further investigate the molecu- 
lar dynamics, an in-depth modeling involving photoab- 
sorption cross sections, Auger decay rates and molecu- 
lar potential surfaces are necessary. Pump-probe exper- 
iments with X-ray pulses much shorter than the Auger 
decay time may in the future provide effective control of 
the involved dynamics, while also serving as a powerful 
FEL beam diagnostics tool. Finally, this work presents a 
quite general approach for treating ionization sequences 
accompanied by break-up of gas phase species, in cases 
of very high photon densities common at FEL sources. 
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Supplemental Material to 'Multiphoton ionization as a "clock" to reveal 
molecular dynamics with intense short X-ray free electron laser pulses' 
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1. Pathways of ionization through photoionization and Auger decays processes 

The photoionization and Auger decays (PA) pathways of N 2 are summarized in the chart 
below (* indicates states with Is hole(s); the bold states and the double arrowed pathways 
are what we are interested in this work; horizontal lines correspond to ionization of one atom 
and vertical lines correspond to ionization of the other atom while the molecule dissociates; 
molecular dissociation starts at the second column): 
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Each PA step removes two electrons. The first PA step mainly doubly ionizes the neutral 
molecules in the ground state. The resultant dicationic ions are in the semi-bound N 2 , 4 " 
states or excited states which lead to dissociation. The latter process dominates following 
the first ionization by the X-ray pulses The asymptotic limit of the dissociation after 
the first core ionization and Auger decay is either N + +N + (the (1,1) channel) or N 2+ +N 
(the (2,0) channel) [2j. Among the observed ions, the N + ions are mainly from the (1,1) 
dissociation [1]. Possible processes leading to the N 2+ are 1) symmetric dissociation of N^ 4 " 
following the PA process on ions in the bound states of N^; 2) symmetric dissociation of 
Ng" 1 " following DCH Auger decay; and 3) (2,0) dissociation or PA process on N from the 
(2,0) dissociation. According to previous experiments 3H5(, process 2) and 3) are minor 
contributions and we consider process 1) as the main cause of the N 2+ ions. Highly charged 
ions are produced by the PA process or photoionization alone on ions of lower charge states, 
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as shown in the chart. For example, during the dissociation of (1,1) channel, the N + ion 
may absorb another photon and is further ionized by the subsequent Auger decay process, 
resulting in N 3+ ; this N 3+ ion can be ionized to higher charge states by further PA processes 
(see columns in the above chart). N 4+ without DCHs and N 5+ cannot ionize further by 
Auger decay process and thus double charge increase by one absorption step is not possible. 
However, sequential photoabsorption can still occur in these fragments as long as they retain 
electrons, in which case the ions N 5+ , N 6+ and N 7+ will be produced. Similarly, the other 
atom in the molecule can be ionized to higher charge states (see rows in the above chart). 

The fast Auger decay mechanism plays an essential role here by refilling the core holes 
in few femtoseconds and allowing for the PAPA. . . sequence to continue. With even higher 
pulse intensities, the core Is orbitals can be depleted of electrons before the Auger decay 
occurs. Such double core hole (DCH) processes reveal novel physics BBS , but they are 
not essential in investigating molecular dynamics and play a minor role in this study. 

2. Extraction of average ion kinetic energy 

To extract the single ion kinetic energy (KE) value we carry out the following procedure. 
First, we model a set of "basis" time-of-flight (TOF) traces by using a well calibrated 
spectrometer geometry in SIMION for each observed charge state of nitrogen atom. Each 
of the TOF traces from the basis set corresponds to a small range of initial KE values of 
the given ion. Next, we take a linear combination of the basis TOF spectra that produces 
the best fit to the experimental TOF spectrum of the corresponding charge state. The 
resulting weights of this linear combination are used to obtain the ion distribution in the 
energy domain. This distribution is further integrated to get the average KE for a given 
charge state. 

3. Model for exacting the timing information of the ionization 

We apply several assumptions to the extraction of the timing information, which do not 
affect qualitatively the results of this work. First, as seen in the chart above, in the column 
of (m,n) with fixed n, the (n,n) state dominates the contribution to n+ ions because they 
require the fewest of photons and the contribution is double of that from the asymmetric 
states in the same column. Due to the latter, even though there is possibility that the (m,n) 
state is populated more than the (n,n) state (along the rows in the chart), the contribution 
from the (n,n) state could be still comparable with other states or even dominates. The 
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above supports the assumption i) in the main text and the validity is confirmed with the 
observations that 1) there is no clear correlation between KEs of different charge states 
which is expected for pure (m,n) contributions where m^n; and 2) for the N + and the N 2+ , 
the measured KEs agree with what is expected for dissociation on the (1,1) and (2,2) curve 
starting at the R e of N 2 (~1.1 A) [8], which is consistent with the scenario of the first PA 
step. Based on the current available data, we cannot extract the kinetic energy distributions, 
from which the exact fragmentation channels may be retrieved. To fully resolve different 
pathways, covariance or coincidence experiments on molecular fragmentation are desired. 
Second, the assumption of time interval between photoionization steps (assumption iv) is 
validated by the X-ray temporal profile 9J. Finally, there are pathways beside the ones used 
that lead to the observed charged ions. We have extracted the time and R using different 
pathways involving extreme asymmetric states of (m,n) with m>(n-3) and find they lead 
to slightly smaller R and time for each PA step compared with results using the assumed 
pathway. The assumptions and approximations applied in our model may add systematic 
errors to the absolute values of time and R, but remain a valid approach for showing the 
relative quantities in the dynamics. 

In the model, the average time between creation of N 4+ and the next photo-absorption 
is Atavgj = (P*p*0"i*^i)~ 1 , where P is the photon flux, p is the target density, cr^ is the 
atomic photoionization cross section, and Ui is the weight of the particular transition. E.g. 
(1,1)— >(3,1) transition would have the weight of 2, because there are two N + ions that can 
be ionized to N 3+ ; while (3,1)— >(3,3) has weight of 1. Given this At avg j, the potential 
corresponding to the given molecular ion, and initial internuclear separation Ri we find the 
increase of the kinetic energy AEj. The summation of all the AEj along the complete multiple 
PA sequence gives the average kinetic energy release (KER) for a given path. The sum of 
KER over all ionization steps gives the average KER for a given charge state. The only 
unknown here is the product of the P*p which is obtained by comparing to the experimental 
data. 

4. Multiple Nl + states 

We only apply the potential curve of 1 A 9 (l7r -2 ) for N 2 + to extract the ionization dynam- 
ics. According to FEL and synchrotron x-ray data, the dominant states that dissociate 
populated by the PA sequence are l A g (lTX' 2 ) (24%), %(l7r~ 2 ) (14%) and 1 n 9 (2a- 1 l7r,- 1 ) 
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(9.8%) 1£)|, QjJ and the motion of ions on these different N| + states may vary. However, 
using these states, respectively, the results do not significantly deviate from that using the 
1 A g state. Therefore, including these states won't change significantly the dynamics map 
and our results show the relative timing and position R of the PA sequences. 
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